double inv_norm(double  y);